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What happens to a neutron star or white dwarf near its maximum mass limit when it is brought 
into a close binary orbit with a companion? Such situation may occur in the progenitors of Type la 
supernovae and in coalescing neutron star binaries. Using an energy variational principle, we show 
that tidal field reduces the central density of the compact object, making it more stable against radial 
\£l • collapse. For a cold white dwarf, the tidal field increases the maximum stable mass only slightly, but 

^\ | can actually lower the maximum central density by as much as 30%. Thus a white dwarf in a close 

. binary may be more susceptible to general relativistic instability than the instability associated with 

electron capture and pycronuclear reaction (depending on the white dwarf composition). We analyse 
the radial stability of neutron star using post-Newtonian approximation with an ideal degenerate 
■ neutron gas equation of state. The tidal stablization effect implies that the neutron star in coalescing 

neutron star-neutron star or neutron star-black hole binaries does not collapse prior to merger or 
tidal disruption. 



PACS Numbers: 97.80.Fk, 04. 25. Dm, 04.40.Dg, 97.60. Jd 



. Introduction 

It is well-established that there exist upper limits to the mass and central density of degenerate compact objects, 
white dwarf (WD) and neutron star (NS) pjM. For WD, the mass asymptotes to the Chandrasekhar limit as the 
central density increases, bounded by general relativistic radial instability, electron capture or pycronuclear reaction. 
For NS, the existence of an absolute mass limit is truly a general relativistic phenomena: pressure, which supports 
, the star against gravity, also acts as a source of gravitation. 

■ WD and NS frequently appear in binary systems. Type la supernova is thought to arise from accreting WD in a close 
binary or the merger of a WD- WD binary || . Of particular interest is the coalescing NS-NS binaries and NS-black hole 
(BH) binaries, which are the most promising sources of gravitational waves that could be detected by interferometers 
such as LIGO and VIRGO @-|). NS binary merger is also considered to be a natural engine that drives cosmological 
gamma-ray bursts. Although the inspiral at larger orbital radius may be treated by post-Newtonian expansion 
technique , coupled with semi-analytic studies of the hydrodynamical effects |sHl3f , quantitative understanding of 
^ . the orbital evolution at small separation and the final merging requires full numerical simulation, which still is in its 
£j ; infancy @||. 

The present study is motivated by the recent general relativistic hydrodynamical simulations of Wilson et al. JnJ 
which reveal evidence that "general relativistic effects may cause otherwise stable neutron stars to individually collapse 
prior to merging" . Obviously, this requires the neutron star to have a mass close to its maximum value to begin with 
(at large separation). However, it does raise a question as to how tidal field modifies the mass limit and central density 
limit of a compact object. In this paper, we analyse the radial stability of neutron star and white dwarf under the 
influence of the tidal field of a companion. We show that the tidal effect generally stablizes the star. We adopt the 
energy variational principle as widely used in the analysis of stability of isolated WD and NS [pclT 



General Consideration 

Let E(p) be the energy of an isolated object of baryon mass M as a function of the central density p. The function 
E(p) consists of internal energy and self-gravitational energy, including their relativistic corrections. The equilibrium 
central density po is obtained from the condition (dE/dp) = 0; radial stability requires (d 2 E/dp 2 ) > 0. 

Now for the star in a binary, the dimensionless tidal distortion is of order ~ e = (M' /M)(R/r) 3 , where R cx (M/p) 1 / 3 
is the mean stellar radius, M' the companion mass, r the orbital separation. The total stellar energy can be written 
as £(p) = E(p) + Wt(p). Note that the tidal distortion modifies the "intrinsic" energy E(p), but we will group this 
correction into Wt(p). Thus the function Wt(p) consists of (i) the correction AW to the self-gravitational potential 
energy, (ii) the interaction energy Wi between M' and the tide-induced quadrupole, and (iii) the kinetic energy T s 
of internal fluid oscillation and rotation. The first two contributions are of the same order ~ (GM 2 / R)e 2 tx 1/r 6 , 
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but have opposite signs, while the third T s ~ MR 2 e 2 ft 2 oc 1/r 9 is a factor (1 + M' /M)(R/r) 3 smaller, and will be 
neglected Q. 

To calculate Wt, we model the star as a polytrope (with equation of state P = Kp 1+1 / n ; n is the polytropic index). 
Approximating the tidaly deformed star as an ellipsoid with axes = + aj) (ai is along the tidal bulge), the 
tidal distortion can be calculated in Newtonian theory. To leading order in e, we have 21 

5 W fR\ 3 5 M' fR\ 3 

where q n — n n (l — n/5), and n n is defined such that the moment of inertia of an isolated object is / = (2/5)n n MR 2 
(k„ specifies the mass concentration within the star). Note that a\ + a 2 + ct 3 = to leading order in e. The two 
contributions to Wt are 

3 GM 2 4 . 2 2 2 

AW = - — — {a 1 + a 2 + a 3 - a x a 2 - a 2 a 3 

5 — n R 45 

. 3 GM' 2 R 5 

-aia 3 ) = -K n qn « , (2) 

and 

GM' ( 2 \ 
W< = -^3- ( -K n MR 2 ) (2a 1 -a 3 -a a ) 



3 GM' 2 R 5 

- --^K n q n -g . (3) 

Thus, up to order e 2 , the total tidal energy is 

Wi = -A^ ~ -A^(M W V3, (4) 

with A = (3/4)k„<2„. Note that the negative sign in Wt is crucial for the tidal stablization of radial mode. 

The equilibrium condition in the presence of tide requires d£ / dp — 0. Thus the density change Sp = p — po (where 
Po is the equilibrium density of the isolated object) due to the tidal field is given by 



Sp = 



5W t (d 2 E 



3p \ dp 



(5) 

Po 



with the expression evaluated at po. Since (d 2 E/dp 2 ) > for stable configuration, we see that 5p < 0, i.e., the tidal 
field reduces the central density of a stable object. If the star is not too close to the stability limit, we have an estimate 
Sp/po ~ (— £ 2 ), which is of second order in the tidal deformation e. Heuristically, a tidally-distorted object is less 
bound gravitationally (compared to a spherical object of the same mass), thus its volume expands in order to satisfy 
hydrostatic equilibrium. 

Now consider how the tidal field changes the maximum mass M m of the object. Let p m o and M m o be the zero- 
tide values of the maximum density and maximum mass, at which dE/dp = d 2 E/dp 2 — is satisfied. Taking the 
difference between (88 / '<9p)| Pm ,M m = and (dE/dp)\ Pm0y M m0 = 0' we obtain 



hW t ( d 2 E V 



SM m = M m -M m0 = ^{—) (6) 



1 We assume the star has negligible viscosity to be tidally synchronized, as is the case for compact objects . We also assume 
the star has zero intrinsic spin — Finite spin tends to stablize the star against radial collapse, and can be treated separately 
(e.g., Q). Note that there is no ambiguity in our definition of the function £(p) for stability analysis: One could add to £(p) a 
term T — T s + T or b associated with the kinetic spin and orbital energies (so that £ can be considered as the total energy of the 
binary system). However, differentiation of T with respect to p under fixed total angular momentum J and fluid circulation C 
gives {dT / dp) j.c = 2T s /3p, and the orbital motion does not affect the stellar structure. 
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bmce we can sho^\(d 2 E/dpdM) < 0, we find 6M m > 0, i.e., the tidal field increases the maximum mass for stability. 
An order of magnitude estimate gives SM m /M m Q ~ Ae 2 . 

Note that two implicit assumptions have been made in our analysis: (i) The star "relaxes" to its equilibrium shape 
even when the density is out of equilibrium. This separation of the radial motion and the "shape adjustment" makes 
our analysis of radial instability more transparent. It is valid because near the stability limit, the radial oscillation 
has almost zero frequency, (ii) The displacement of a fluid element inside the star is a linear function of the fluid 
position. This is exact only for the n = 3 polytrope and in the incompressible limit, but otherwise corresponds to an 
approximate trial wavefunction in the variational principle. 

White Dwarf 

We now consider the stability of white dwarfs. Near the maximum mass, the stellar density profile resembles that 
of a n — 3 polytrope. The "intrinsic" energy can be written as jl],|| 

E = E mt + W + AE int + AE GR , (7) 

where Ei nt ~ Mp 1 ^ 3 is the internal energy of the ultra-relativistic electrons, AEi nt ~ Afp -1 / 3 is the cor- 
rection due to finite electron mass, W ~ —M 2 /R ~ — M 5 / 3 p 1 / 3 is the potential energy of self-gravity, and 
AEgr ~ — (M 2 /R)(M/R) ~ —M 7 1 3 p 2 / 3 is the general relativistic (post-Newtonian) correction. For isolated cold 
WD, this gives for the maximum central density and maximum mass p m o = 2.737 x 10 10 gem -3 , M m o = 1.4156M , 
the correcponding minimum radius R m o — 1110 km (We assume that the number of electron per nucleoli is Y e = 0.5). 

We parametrize the strength of the tidal field by the dimensionless ratio (3 = (M' /M m o)(R m o/r) . To avoid tidal 
disruption we require /3 <^ 0.1. With A ~ 0.01 (for n = 3 polytrope) we then have 

Expressing mass in the units of Mq, density in 10 10 gem -3 , energy in 10 51 erg, the total energy can be written as 

£ = 3.7129Mp 1/3 - 2.8895M 5/ V /3 + 0.0457Mp -1/3 

-0.0105Af 7/ V /3 - 0.14 f3 2 (M/p) 5/3 . (9) 

Solving (d£/dp) — (d 2 £ /dp 2 ) = yields the maximum density and maximum mass as a function of j3. The result 
is shown in Figure 1. The increase in the maximum mass is rather small (^ 0.1%). However, the maximum density 
decrease can be rather substantial, reaching as much as ~ 30%. For small 0, we have SM m /M m Q ~ 0.05/3 2 , and 
Spm/pmo ~ — 17/3 2 . Figure 1 also shows the central density of three constant-mass sequences with mass slightly 
below M m o. We see that as the binary separation decreases ([3 increases), the central densities also decrease, and 
always remain smaller than the maximum density p m {P) allowed for radial stability. 

The fact that the maximum density p m {P) for stability decreases with decreasing r may have some interesting 
astrophysical consequences. The maximum central density of a normal Carbon-Oxygen white dwarfs is set by the 
threshold of electron capture on 16 0, at p cap — 1.9 x 10 10 gem -3 (the pycronuclear reaction between 12 C may set 
in at smaller density, of order 10 10 gem -3 ; but this increases as the 12 C abundance decreases ^2|). From Fig. 1 we 
see that p m {f3) can drop below p ca p, therefore GR effect becomes more important than neutronization. This implies 
that an accreting white dwarf in close binary can be more susceptible to collapsing to NS by general relativistic radial 
instability. Of course, the real situation may be complicated by the finite temperature due to the high accretion rate. 

Neutron Star 

The equation of state of nuclear matter is uncertain. For the purpose of illustration, we adopt EOS to be that of 
an ideal degenerate neutron gas and analyse the radial stability in the post-Newtonian approximation. Assuming the 
stellar density profile to be that of the n — 3/2 polytrope, the "intrinsic" energy can still be written as in Eq. (7), 
except that Ei nt ~ Mp 2 / 3 is the internal energy of nonrelativistic neutron, AEi nt ~ —Mp 4 / 3 is the correction due to 
special relativistic effect. The maximum baryon mass of an isolated NS bassed on this model is M m o — 1.11O8M , 



2 Consider a sequence of stellar models parametrized by the central density. Differentiating dE/dp = we obtain: (d 2 E /dp 2 ) + 
(d 2 E/dpdM)(dM/dp) = 0. The stable branch of the sequence satisfies (d 2 E/dp 2 ) > and dM/dp > 0, while the unstable 
branch satisfies (d 2 E/dp 2 ) < and dM/dp < 0. In both cases we have (d 2 E/dpdM) < 0. 
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the maximum central density p m Q — 7.415 x 10 15 gem -3 , and the corresponding minimum radius R m a — 10.48 km. 
We again parametrize the tidal strength by the ratio (3 = (M' /M m o)(R m o/r) . The tidal energy can then be written 
in the form of Eq. (8) except that we use A ~ 0.1 appropriate for n = 3/2 polytrope (The uncertainty can be readily 
absorbed into the definition of (3). The GR correction to the tidal potential energy is of order (M'/r) smaller, and is 
neglected. Expressing mass in the units of M , density in 10 15 gem -3 , and energy 10 53 erg, the total energy function 
can be written as 

£ = 0.85168Afp 2/3 - 1.5968M 5/ V /3 - 0.02887M/? 4 / 3 

-0.16774M 7/3 p 2/3 - 7.356 (3 2 (M/p) 5/3 . (10) 

Figure 2 shows the maximum density p m and maximum mass M m as a function of (3. We see that, as in the white 
dwarf case, p m decreases and M m increases with increasing f3. For small /3 (^ 0.15), we have 

SM -- 15W < 0.33 (11) 



M m0 ~ 4(2E m + W) 

1 RWl 

= -2.7/3 2 . (12) 



S Pm _ 5W t {28E in + 15W) 



PmO 

2W(2E m + W) 

Figure 2 also depicts the densities of constant-mass sequences with M slightly below M m o. Again, we see that these 
densities p(f3) always remain below p m (/3), i.e., these stars are stable against radial perturbation until the tidal limit 
(3 ~ 0.23 is reached. 

Another way to look at the tidal effect on the radial stability is to construct sequence of stellar models with varying 
central densities at a given fixed (3. Figure 3 shows the M — p curves of several such sequences. Only models in the 
dM/dp > branch are stable, and the maximum density p m and maximum mass M m are given by the values at the 
turning point dM/dp = 0. One can readily see that the effects of the tidal field is to raise M m and lower p m . For 
[3 <; 0.23, no stable configuration exists, and this corresponds to the tidal disruption limit. 

Discussions 

The above analysis demonstrates that tidal effect tends to stabilize neutron star and white dwarf against radial 
collapse, at least within the framework of post-Newtonian theory. Our method is accurate for white dwarfs, but only 
approximate for neutron stars. However, we think it is unlikely that higher order GR corrections or the use of more 
sophisticated nuclear equation of state |19| will change our qualitative results for neutron stars, although the precise 
numbers can certainly change. 

In the cases of NS-NS binaries, one might still consider the possibility of neutron star collapse prior to merging 
when there is stable mass transfer from its lower mass companion. However, by the time mass transfer starts, the 
binary must already have encountered the orbital dynamical instability as a result of strong tidal interaction [fl0|| . 
This tidal instability, enhanced by general relativistic effects (2^] leads to rapid coalescence of the binary within a few 
orbits HQ. 

I thank Kip Thorne for comments on an earlier version of the paper and helpful suggestions. This research is 
supported by the Richard C. Tolman Fellowship in theoretical astrophysics at Caltech, and NSF Grant AST-9417371 
and NASA Grant NAG 5-2756. 
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FIG. 1. The maximum central density (heavy solid line) and maximum mass (dashed line) of white dwarf as a function of 
the tidal ratio f3 = (Af'/M m o)(-Rmo/r) 3 . The solid lines show the central densities of three sequences with mass slightly smaller 
than maximum mass of a isolated WD M m o = 1.4156-Mq. 

FIG. 2. The maximum central density (heavy solid line) and maximum mass (dashed line) of neutron star as a function of 
the tidal ratio /3 = (M 1 /M m o)(Rmo/r) 3 . The solid lines show the central densities of three sequences with mass slightly smaller 
than maximum mass of an isolated NS M m o = I.IIO8M0. 

FIG. 3. The mass-central density curves of sequences of NS models with constant values of /3. (j3 = corresponds to isolated 
NS). 
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